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VUV  FOURIER  TRANSFORM  SPECTROSCOPY  OF  THE  (0,0)  BAND  OF  NO 


« 


K.  Yoshino,  J.R.  Esmond,  J.E.  Murray  A.P.  Thome 

Y.  Sun,  A.  Dalgarno,  W.H.  Parkinson 


Harvard-Smithsonian  Center  for 
Astrophysics 
Cambridge,  MA  02138 


Blackett  Laboratory 
Imperial  College 
London  SW7  2BZ,  UK 


Using  the  VUV  Fourier  transform  spectrometer  (FTS)  at  Imperial  College  we  have  observed 
the  (0,0)  and  (1,0)  bands  of  NO,  at  191  nm  and  183  nm  respectively,  in  absorption.  A  high 
current  hydrogen  lamp  was  used  as  the  background  continuum  source  with  a  0.3  m  grating 
spectrometer  as  a  bandpass  filter.  The  absorption  spectrum  was  observed  at  various  pressures 
of  NO  ranging  from  0.004  to  0.089  torr,  and  temperatures  of  295  K  and  78  K.  To  fully  resolve 
the  sharp  NO  lines  (Doppler  width  0.11  cm  *  at  295  K)  we  chose  an  FTS  resolution  of  0.06  cm'* 
(2  mA).  We  report  here  on  the  first,  absolute,  cross  section  measurements  made  for  the  (0,0) 
band  along  with  new  improved  energy  precision  for  the  low  lying  rotational  levels  (up  to 
1=201/2  for  the  C^II  and  1=81/2  for  the  B^II  electronic  states).  A  deperturbation  procedure  to 
analyse  the  energy  level  structure  will  be  presented. 

This  work  is  supported  by  NSF  Division  of  Atmospheric  Sciences  grants  ATM-91-16552  and 
ATM-90-19188  to  Harvard  College  Observatory. 
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VUV  Fourier  Transform  Spectroscopy  of  the  ^(0,0)  Band  of  NO 
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Photo-destruction  of  NO  is  an  important 
stratospheric  process 

175  -  205  nm  penetration  controlled  by  O2 

Atmospheric  model  computations  involving 
O2  and  NO  must  be  done  on  a  line  by 
line  basis 

Such  measurements  have  been  completed  for  most 
of  the  Schumann- Runge  band 

To  fully  resolve  the  linewidths  of  the  NO  band 
we  require  resolving  powers  approaching  10® 


8(0,0)  AND  PiT.O)  BANDS  OF  NO 
AND  (5.0)  S-R  BAND  OF  Oo 
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Czerny- Turner  Spectrometer 
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Wavenumber  and  term  differences  between  this 

work  and  previous  studies 
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Table  1 


Determination  of  the  mixed  C  and  B  Tye  terms 
using  rotational  lines  from  different  branches 


J 

Rll 

Pll 

R12 

P12 

Qll 

Q12 

Mean 

std 

0.5 

52372.72 

52372.69 

52372.70 

0.02 

1.5 

52380.29 

52380.31 

52380.32 

52380.29 

52380.30 

0.01 

2.5 

52391.00 

52390.99 

52391.00 

52390.99 

52390.95 

52391.02 

52390.99 

0.02 

3.5 

52402.34 

52402.33 

52402.27 

52402.31 

0.04 

4.5 

52414.02 

52414.04 

52414.03 

0.02 

5.5 

52430.65 

52430.65 

52430.63 

52430.64 

52430.64 

0.01 

6.5 

52453.45 

52453.45 

52453.45 

52453.44 

52453.45 

0.00 

7.5 

52480.79 

52480.80 

52480.76 

52480.78 

0.02 

8.5 

52512.26 

52512.26 

52512.25 

52512.25 

52512.26 

0.00 

9.5 

52547.75 

52547.70 

52547.76 

52547.74 

0.03 

10.5 

52587.19 

52587.19 

52587.24 

52587.26 

52587.22 

0.03 

11.5 

52630.65 

52630.63 

52630.64 

0.01 

12.5 

52678.06 

52678.06 

52678.04 

52678.06 

0.01 

13.5 

52729.40 

52729.39 

52729.40 

0.01 

14.5 

52784.74 

52784.72 

52784.74 

52784.73 

0.02 

15.5 

52843.98 

52844.02 

52844.00 

0.02 

16.5 

52907.15 

52907.16 

52907.16 

0.00 

17.5 

52974.30 

52974.35 

52974.33 

0.03 

18.5 

53045.38 

53045.39 

53045.38 

0.01 

19.5 

53120.39 

53120.39 
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The  eigenstates  of  the  NO  molecule  are  mixtures  of  four  Hund’s  case  (a) 
basis:  C'^ni/2(t'  =  0,  J),  =  0,  J),  -  7,  J),  and  B^Tlz/2[v  = 

7,  J).  The  corresponding  4x4  Hamiltonian  is  given  in  Table  2a.  Least  square 
nonhnear  fitting  of  the  observed  enenrgy  values  to  the  Hamiltonian  results  in 
the  molecular  constants  (Table  2b).  In  the  fitting,  we  fixed  the  D  constant 
for  the  B^H  at  4.68  x  10~®  (calculated  using  the  RKR  potential  of  Gulluser 
and  Dressier). 

Table  2a:  Hamiltonian  for  each  block 


!  ’n.rt 

^  ^n3/2 

■  “Hi/o  I  T,  -  .4/2  +  B{X  -f  1)  -  D\{X  -r  1)^  ^ 

X.  symmetric 

-tv{l  -  +  q{X  ~2-  2v/r' 

-l)/2' 

!  ^Hs/j  i  ^X[-B  -r  2DX 

-  .4/2  B{X  -  1  j 

j  i  -p/4  ^  qi±VX  -rl-l)/2] 

X  =  (J  —  1/2)(J  4-  3/2),  upper  (lower)  signs  are  for  e(f)  parities. 

{Bma\H\Cmu')  =  HbcSuw 


Table  2b:  Molecular  constants  in  cm  units 


■  C2n(u  =  0) 

T 

.4 

B 

D 

P 

q 

this  work 

52371.26 

3.21 

1  1.9945 

6.2  xlO-'^  : 

-0.0059  : 

-O.OIG 

A&V 

1 _ 

3.1957 

i  1.994555 

5.8694x10-“  ! 

-0.0057  i 

-0.0159484 

B2n(t;  =  7) 

T. 

A  i 

B  ! 

D  1  p 

q 

this  work 

1 

52367.14 

42.31  i 

- . 

1.0320  1  4.68x10-'^  1  0.024 

I  -0.0034 

Hbc  =  5.29,  rms=2.1xl0 

The  energy  origin  is  the  lowest  X^H  level. 
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J(J+1) 


Band  oscillator  strengths  of 
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Atmospheric  Transmission  Issues 
for  the  Air  Force  Airborne  Laser  (ABL) 
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W.  J.  Schafer  Associates  Inc 
1901  N.  Fort  Myer  Dr. 
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(703)  558-7900 


Atmospheric  Effects  on  ABL  Performance 
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Future  work 


Purpose  of  This  Briefing 
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AIRBORNE  LASER  LABORATORY 
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ABL  Propagation  Geometry 
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Laser  T  racker/Beacon  1 .06  1 .3 
High  Energy  Laser  1 .32  others 


ABL  RANGE  vs.  THREAT  CASE 


ABL  LETHAL  ENVELOPE  EXPANDS  WITH  INCREASING  RANGE  OF  TBM  THREAT 

(CS60  =  B470  s  A350) 


ABL  Cases  A,  B,  C,  D 
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Tilt-Only  Compensation  Higher  Turbulence 

Recent  Volcanic  Aerosols  Near-term  Systems 

SAGE  Satellite  Data  Others 


Atmospheric  Transmission  is  One  Loss  Factor 
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INPUT  ASSUMPTIONS  DEVICES/AO 
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Atmospheric  Effects  on  ABL  Performance 
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Future  work 


Total  Atmospheric  Transmission 

Is  The  Product  of  Four  Factors 
COIL  (1.32  n)  Case-A  Path  Example  (12.5/15.3/320) 
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Low-Resolution  (LOWTRAN)  -  Case-) 

Bkg  Strat/Mod  Vole  MLS  12.5/15.3/320  km 
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High-Resolution  Atmospheric  Absorption 
Spectra  near  the  Iodine  Laser  Frequency 
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Wavenumber  (1/cm) 


Sources  of  Aerosol  Extinction 

(Above  Clouds) 
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Satellite  Limb  View 


I/I/JV/L  RECENT  VOLCANOES 
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iMI iCJL  VOLCANIC  AEROSOL  TR ANSMITT ANCE 

CASE  A 


*  Absorption  from  these  Aerosols  could  exacerbate  thermal  blooming. 


A  Brief  Description  of  the  Life  of  a  (Volcanic) 

Stratospheric  Aerosol 
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Atmospheric  Effects  on  ABL  Performance 
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Future  work 


Profile  of  the  Scattering  Ratio  Taken  at  the 
Maximum  -  About  8  Months  After  Pinatubo 
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H.  Jager,  The  Pinatubo  Eruption  in  Relation  to  the  El  Chichon  Event  and  the  Stratospheric 

Background  Load.  0SA  Topical  Conference  on  Remote  Sensing  of  the  Atmosphere,  March,  1 993.  wAT- 


NASA  LIDAR  Backscatter  Results 

1/e  recovery  time  of  about  7.3  months  after  volcanic  aerosol  peak 

Peak  is  reached  about  8  months  after  eruption 

Extinction  varies  with  geographic  location  and  specific  volcano 


M.  Osborn  et  al,  Evolution  of  the  Pinatubo  Volcanic  Cloud  Over  Hampton.  VA.  OSA  Topical 
Conference  on  Remote  Sensing  of  the  Atmosphere,  March,  1993. 


NASA  LIDAR  Backscatter  Results 

1/e  recovery  time  of  about  7.3  months  after  volcanic  aerosol  peak 
Peak  is  reached  about  8  months  after  eruption 


NASA  LIDAR  Backscatter  Results 

1/e  recovery  lime  of  about  7.3  months  f  ft®!  aerosol  peak 
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Atmospheric  Effects  on  ABL  Performance 
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Future  work 


SAGE  SOLAR  OCCULTATION 
MEASUREMENT  GEOMETRY 


1 1 


Latitudinal  coverage  from  79°S  to  79°N.  Provides  15  sunrise  and  15  sunset 
measurements  each  day  whose  latitudes  varied  daily  providing  near  global 
coverage  in  3  -  4  weeks. 


SAGE  (STRATOSPHERIC  AEROSOL 
AND  GAS  EXPERIMENT) 
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Orbit:  Tailored  to  maximize  geographic  coverage  of  solar 

occultation  measurements. 
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SAGE  1990  -  91  Middle  East 


SAGE  1984  Middle  East 
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Aerosol  Extincllon  Coefficients  (km 
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SAGE  1986  Middle  East 


426 


Aerosol  Extinction  Coefficients  (km 


SAGE  1987  Middle  East 
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Aerosol  Extinction  Coefficients  (krn 
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SAGE  1991  Middle  East 
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Atmospheric  Effects  on  ABL  Performance 
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Future  work 
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ABL  Performance  Envelope  Using  LIDAR  & 
Satellite  Measured  Aerosol  Extinction  Profiles 
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100  200  300  400  500  600  700  800 

Horizontal  Range  (km) 
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Impact  of  Pinatubo  Aerosols;  Bl  Scenario 
1  Baseline  ABL  in  Loiter  Orbit  - 132  targets 
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Atmospheric  Effects  on  ABL  Performance 
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Future  work 


Volcanic  Aerosols  -  Issues 
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summer)  will  show  more  detail  on  recovery  rates. 

Impact  of  backscatter  on  laser  tracker  and  A/0  beacon 
needs  to  be  assessed. 


Wednesday  9  June  1993  a.m. 


SESSION  D:  Non-LTE  SPECTROSCOPY  APPLICATIONS 
Chair:  Richard  H.  Picard,  PL/GPOS 


440 


SHARC 

A  MODEL  FOR  CALCULATING  ATMOSPHERIC  RADIATION  UNDER 
NON-EQUILIBRIUM  CONDITIONS 


David  Robertson,^  Robert  Sundberg,*  James  Duff/ 

John  Gruninger/  Steve  Adler-Golden/  Ramesh  Shanna/ 
and  Rebecca  Healey^ 

^Spectral  Sciences,  Inc.  ^Phillips  Laboratory/GPOS 

99  South  Bedford  Street,  #7  29  Randolph  Road 

Burlington,  MA  01803-5169  Hanscom  AFB,  MA  01731-3010 

^Yap  Analytics 
594  Merrett  Road 
Lexington,  MA  02173 


SHARC  was  developed  by  the  Air  Force  to  provide  both  research  and  systems-level 
predictions  for  atmospheric  IR  radiance  for  arbitrary  paths  within  the  50  to  300  km  altitude 
regime  and  in  the  2-30  ptm  spectral  region.  The  code  calculates  LTE  and  NLTE  emissions  from 
the  significant  atmospheric  IR  radiators,  CO^,  NO,  O3,  H2O,  CO,  OH,  and  CH*.  Molecular 
excited  state  populations  are  calculated  with  a  Monte  Carlo  model  for  layer-layer  radiative 
excitation  and  energy  transfer  with  a  flexible  chemical  kinetics  module  derived  from  the 
CHEMKIN  Code  developed  by  Sandia,  Livermore.  Radiation  transport  calculations  are  based 
on  an  equivalent-width  line-by-line  (LBL)  approach  with  a  spectral  resolution  of  about  0.5  cm'^ 
The  LBL  algorithm  uses  the  Roger-Williams  and  Curtis-Godson  approximations  for  the 
equivalent  widths  of  combined  Doppler-Lorentz  (Voigt)  lineshapes.  Spectroscopic  data  are  taken 
from  the  HTTRAN  line  atlas  which  has  been  augmented  with  additional  O3,  NO  and  NO"^  lines. 
SHARC  also  has  an  auroral  module  that  describes  electron  dosing  and  solves  the  time/energy 
dependent  rate  equations  to  calculate  secondary  electron  distributions  and  the  resulting  IR 
emissions  from  CO2,  NO  and  NO"^.  This  module  is  fully  embedded  in  the  ambient  part  of  the 
code  so  that  radiance  calculations  for  paths  passing  through  a  finite  auroral  region  are  possible. 
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SHARC 

A  MODEL  FOR  CALCULATING  ATMOSPHERIC  RADIATION 


Lu 

Q 

CO 

Uj 


CO 

cr 

Uj 

CD 

I 

CO 


X3 


q: 

LU 

CQ 

O 

cc 


o 

CO 

cr 

Lu 

cn 

CD 

cr 


Lu 

CD 

-u 

Q 

CO  • 

1  CO 

cr  ^ 

Lu 

-U 

Q  . 
-cc  CO 
Lu 
Lu  CO 
^ 
Lu  Lu 

CO  Co 
CO 


cr 

Uj 

CO 


CD 

cr 

CO 


cr 

h- 

CO 

Lu 

cr 

CO 


3: 

CD 

r) 


CD 


CO 

.:r 

>-i 

Lu 


-J  . 
CO 


Lu  Co 

r:  h-i 


<c 

CO 

Co 

Lu 

cr 

Lu 

cr 


>- 

-j 


'<5: 


CO 


to 


Lu 

CO 

Lu 

cr 

Lu 

CD 

CO 


c_ 

Q. 


Lu 

CD 

CD 

5: 


CO 

CO 


T3 

Qj  Lu 

■40  t~-1 

C 

QJ  Lu 
cn  a; 

03  ^ 


§ 

5 


-u 

CD 


Lu 


cr 

Lu 

r 

cr 

CO 

CD 


442 


9  JUNE  1993 


ACKNOlf^LEDGEMENTS 


• — I 

a 


2: 

in 

Q 

o 

^-1 

Uj 

h- 

CD 

1 

1-1 

Nj 

U. 

^-l 

:s: 

u> 

a. 

CD 

<3: 

CD 

Uj 

QC 

cT 

<D 

Q 

h- 

Uj 

cn 

CO 

>- 

h- 

Uj 

CD 

t— 1 

U. 

Uj 

Q 

CD 

Uj 

Lu 

q: 

Uj 

Q 

CO 

cn 

Uj 

CO 

Q 

u. 

CO 

CL 

Uj 

►-1 

cn 

CD 

cn 

CD 

CD 

1— 1 

t-i 

CD 

h~ 

Uj 

CO 

h~ 

1-1 

cc 

Q 

CD 

J— . 

'<5: 

CO 

CO 

CD 

Uj 

Lu 

3: 

d: 

<D 

QC 

h- 

1— . 

CO 

Uj 

CO 

o 


QC 

CD 

CO 

CO 

Lu 

CJ 

Uj 

CD 

Uj 

QC 

Q. 


443 
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SHARC  OBJECTIVES 
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MODEL  ATMOSPHERIC  EMISSION  PROCESSES 

SUPPORT  AF  MEASUREMENTS  PROGRAMS 
CIRRIS  lA.  SPIRIT  II.  MSX 

VALIDATION  ^ITH  FIELD  AND  LAB  DATA 


ILLUSTRATIVE  SHARC  CALCULATION 
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AURORAL  MODULE  \tilTH:  NO-h,  NO.  CO2 

EQUIVALENT -yilDTH  LBL  RADIANCE  CALCULATIONS 
MULTIPLE  REGIONS  (CHANGING  ATMOSPHERIC  CONDITIONS) 
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RADIATIVE  EXCITATION  MODULE 

(NEMESIS) 
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SHARC  AURORAL  MODULE 
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ILLUSTRATIVE  AURORAL  ENHANCEMENT 
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Frequency  (cm 


COMPARISON  TO  AURORAL  FIELD  DATA 
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STRONGLY  PREDOSED  Frequency  (cn 


SPECTRAL  RADIANCE  MODULE 
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50-70  TIMES  FASTER  THAN  TRADITIONAL  LBL  APPROACH 


SPIRE:  DAYTIME  CO^  (4.3  liin) 
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SHAFIC-3  TERMINATOR 


456 


In-Band  Radiance  (W/cm^/sn) 
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Limb  Radiance  (W/cm^/sr) 


FUTURE  PLANS 
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SPECTRAL  REGIONS 


SUMMARY 
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AVAILABLE  FROM  PL/GPOS  -  CONTACT 

-  DR.  RAMESH  SHARMA  (617)  377-4198 


NON-LTE  STUDIES  OF  THE  15-^m  BANDS  OF  CO^  FOR 
ATMOSPHERIC  REMOTE  SENSING 


David  P.  Edwards/  Manuel  Ldpez-Puertas/  and 
Miguel  Angel  L6pez-Valverde^ 


^National  Center  for 
Atmospheric  Research 
Boulder,  Colorado,  USA 


^Instituto  de  Astrofisica  de 
Andalucia 
Granada,  Spain 


The  new  line-by-hne  non-LTE  calculation  capability  of  the  GENLN2  radiative  transfer  code  is 
described.  Non-LTE  model  implementation  and  molecular  state  vibrational  temperature  input 
requirements  are  discussed  for  studies  of  the  15-pm  bands  of  CO^.  Monochromatic  and  band- 
integrated  radiance  calculations  have  been  performed  for  atmosphenc  limb  view  tangent  heights 
between  50  and  120  km  for  non-LTE  night  and  daytime  conditions.  Non-LTE  radiance 
considerations  are  shown  to  be  important  for  the  15-pm  CO^  bands  for  tangent  heights  greater 
than  70  km,  the  magnitude  of  the  divergence  from  LTE  values  and  diurnal  variation  bemg 
dependent  on  the  band  and  kinetic  temperature  profile.  We  show  the  importance  of  including 
Lorenzian  line  wings  and  overlapping  lines.  Calculations  of  synthetic  radiance  spectra  are 
presented  showing  the  non-LTE  effect  for  two  CO^  temperature  sounding  channels  of  inst^ments 
aboard  the  Upper  Atmosphere  Research  Satellite  as  a  demonstration  of  the  model  capability. 
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CRITICAL  TESTS  OF  NON-LTE  RADIATIVE  MODELS  AGAINST 
HIGH-LATITUDE  ROCKET  DATA 


-V 


R.H.  Picard/  J.R.  Winick/  U.  Makhlouf/  AJ.  Paboojian/ 

AJ.  Ratkowski/  K.U.  Grossmann/  D.  Homann/  and  J.C.  Ulwick' 


^Phillips  Laboratory/  Geophysics 
Optical  Environment  Division 
Hanscom  AFB,  MA  01731 


^Stewan  Radiance  Laboratory 
Utah  State  University 
129  Great  Road 
Bedford,  MA  01730 


^ARCON  Corporation 
Waltham,  MA  02154 


%'nivcrsity  of  Wuppertal 
Wuppertal,  Germany 


We  have  earned  out  tests  of  non-LTE  radiation  models  against  field  data  under  conditions  in 
which  the  state  of  the  atmosphere  was  very  well  characterized  and  subject  to  extreme  differences. 
The  data  were  obtained  during  a  series  of  rocket  launches  from  northern  Scandinavia,  supported 
by  ground-based  observations.  Tnese  included,  the  series  of  SISSI  flights,  the  first  of  which  (6 
Mar  1990)  was  part  of  the  DYANA  Campaign,  and  the  MI-1  rocket  launched  10  Feb  1984  during 
the  MAP/WINE  Campaign.  The  measured  state  of  the  atmosphere,  including  temperature  and 
atomic-oxygen  profiles,  is  input  to  the  line-by-line  ARC  (Atmospheric  Radiance  Code)  non-LTE 
code  and  associated  models,  and  model  predictions  are  compared  with  spectral  data  from  a 
rocketbome  Eben-Fastie  spectrometer,  emphasizing  the  CO,  15um  data  We  conclude  that  the 
models  are  able  to  predict  non-LTE  spectral  radiance  very  well  when  provided  good  input  data 
on  the  atmospheric  state.  We  also  show  that  very  significant  differences  between  model 
predictions  and  point  measurements  can  occur  when  climatologies  are  used  to  generate  inputs  for 
the  radiative  models. 
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Critical  Tests  of  Non-LTE  Radiative  Models 
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SISSI  PROGRAM 
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Solar  Depression  Angle  near  5°  all  flights 
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Vertical  probe  -  Comparison  to  limb  scans 
Well  characterized  atmosphere;  crucial  for  testing  IR 
model  response  to  atmospheric  variability 


SISSI-1  (DYANA  Campaign) 

ESRANGE,  Sweden  (68°N,  21  °E) 

6  Mar  90,  0440  UT 

Rocket  -  Skylark  VI 

Apogee;  200km 

SZA:  =95“ 


M-11  (MAP/WINE  Campaign) 

ESRANGE,  Sweden 
1 0  Feb  84,  041 2  UT 
Rocket  -  Skylark  VI 
Apogee:  179km 
SZA:  106“ 
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1  » 


SPLIT  NOSE  CONE 


INFRARED  SPECTROMETER 
(Univ.  of  Wuppertal) 


NITRIC  OXIDE  SENSOR 
(USU/PL) 

FARADAY  ROTATION 
cruG) 

IR  RADIOMETER 
(USUfPL) 

ATOMIC  OXYGEN  SENSOR 

(USUfPL) 

PLASMA  FREQ-  PROBE 
(USUfPL) 

-  BATTERIES 


ATTITUDE  CONTROL  SYSTEM 


RECOVERY  UNIT 
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ARC  (Atmospheric  Radiance  Code)  [cont.] 

3  modules:  .  , 
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CO2  15-Mm  RADIANCE  MODEL  (ARC) 
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Atomic  oxygen  excitation  /  quenching  especially 
important  because  of  large  value  of  rate  constant: 
ko=6.1x10'^^  cm^/s  at  T=300K,  varies  as  (Sharma  & 
Wintersteiner,  1990;  Wintersteiner  et  al.,  1992) 
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ARC  Model  Vibralioiiarrcinperntures  -  SlSSl- 1 
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01-01 


Wavelength  Data  ARC  Model 


ance  ( W/ciii2/.si7iiiti) 


Self-consistent  One-dimensional  Photochemical  Model 

U.  Makhlouf\  J.  R.  Winick^  and  R.H.  Picard^ 


.9-  CO 

b  2 

CO  3 
I— I 


0  0 
C  "D 
0  O 


CD  s 

•a  : 
o  • 

E  X 

od 

o  -CO 
Q.O 

l" 

C 


^  c/5 
<35  Sfi 


C  05 
CD  -C 
■a  jD 

o  -a 

g-s 

T3  O 

go 

°>s 

c  _2 

^  C/) 

o  T3 
C  CCS 

s 

"55  <0 

■a  CO 
O  CD 

C  ZJ 

g 

S  2 

o  *- 

Q.  <D 
O  «5 

t§ 

CCS 

c  O 

O  CCS 
CCS  ‘- 

^  -a 

~  CD 
>  TS 


X  ^ 
O  CD 


0  E 

O)  ^ 
§  § 
O  ° 


0  ^ 

0  2 
O  CD 

£  E 

0  ® 


CO 

E  p 


‘o)  ^ 

o  -Q 
CCS  2 

.i  55 
cc  2 

kx.  k,_ 

z  CO 

E  :i 
^  > 

CD  ^ 


0  ^ 

0  = 
■ii-  O 

0  <^ 
»-  _l 

1  2 
0  CO 
«  CD 

C  ^ 

.9  E 

H— »  O 

^  -O 
0  o 

0  X 

.9  :§■ 


486 


Stewart  Radiance  Laboratory,  Bedford  MA  01730 
Phillips  Laboratory/GPOS,  Hanscom  AFB,  MA  01731 


Atomic  Oxygen  Density;  SISSI  FI 
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Atomic  Oxygen  Density  (cm-3) 


(HJO  sonetpe^i 
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Solar  Zenith  angle  (degrees) 
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ATMOSPHERIC  STRUCTURE  SIMULATION:  AN  AUTOREGRESSIVE 
MODEL  FOR  SMOOTH  GEOPHYSICAL  POWER  SPECTRA  WITH 
KNOWN  AUTOCORRELATION  FUNCTION 


James  H.  Brown 


Phillips  Laboratory/  Geophysics 

29  Randolph  Road 

Hanscom  AFB,  MA  01731-3010 


Within  a  defined  domain,  geophysical  phenomena  often  are  characterized  by  smooth  continuous 
power  spectral  densities  having  a  negative  power  law  slope.  Frequently,  Fourier  transform 
analysis  has  been  employed  to  generate  synthetic  scenes  from  pseudorandom  arrays  by  passing 
the  stochastic  data  through  a  Fourier  filter  having  a  desired  correlation  structure  and  power 
spectral  dependency.  This  paper  examines  the  possibility  of  producing  synthetic  structure  by 
invoking  autoregression  analysis  as  contrasted  with  the  Fourier  method.  Since  computations  that 
apply  multidimensional  fast  Fourier  transforms  to  large  data  arrays  consume  enormous  resources 
and  time,  the  goal  of  this  study  is  to  seek  an  alternative  method  to  reduce  the  computational 
burden.  Future  editions  of  the  Phillips  Laboratory  Strategic  High  Altitude  Atmospheric 
Radiance  Code  (SHARC)  will  feature  an  ability  to  calculate  structured  radiance.  The  methods 
explored  herein  provide  a  process  that  can  compliment  or  in  some  cases  supplement  methods 
presently  being  used. 
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Figure  Captions 


1 .  Expression  for  the  linear  difference  equation.  Future  values  of  a  time 
series  predicted  from  past  values  and  random  data.  So-called  ARMA 
model. 

2.  Expression  for  power  spectra  given  an  ARMA  linear  predictor  model. 
The  a's  and  b's  are  the  AR  and  MA  predictor  coefficients. 

3.  A  theoretical  model  for  power  spectral  density  and  autocorrelation 
function  for  constant  power  law  slope. 

4.  Expression  for  the  coherence  length  parameter  .  Function  of  power 
law  slope  and  parameter  "a". 

5.  Sample  model  PSD  and  autocorrelation  function.  Slope  =-3..  = 

0.2,  Lc  =  5.0  and  a  =  0.05.  Upper  left  quadrant  log-log  PSD  vs  spatial 
frequency,  upper  right  quadrant  linear  PSD  plot,  lower  left  quadrant 
autocorrelation  function  vs  lag,  lower  right  quadrant  power  curve. 

6.  White  noise.  Sample  Gaussian  random  number  sequence,  mean  =  0,, 
standard  deviation  =  0.055,  spacing  100m. 

7.  Histogram  of  sample  Gaussian  random  number  sequence.  Theoreti¬ 
cal  mean  =  0.,  theoretical  S.D.  =  0.055. 

8.  Left  panel,  theoretical  PSD  input  (unmarked),  autoregressive  predic¬ 
tor  PSD  (marked  by  x),  and  simulated  PSD  (marked  by  small  square). 
Right  panel,  corresponding  autocorrelation  functions.  Lq  =  1.75  Km, 
S  =  -5/3,  a2  =  1.02E-03,  spacing  =  100m.  Theoretical  autocorrelation 
function  modified  at  lag  =  0.  Six  predictor  coefficients. 

9.  Autoregressive  simulated  structure  sequence.  Lq=  1.75  Km,  S  = 

-5/3,  a2  =  1.02E-03. 

10.  Histogram  of  autoregressive  simulated  structure  sequence.  Lq  =  1.75 
Km,  S  =  -5/3,  a2  =  1.02E-03. 

11.  Left  panel,  theoretical  PSD  input  (unmarked),  autoregressive  predic¬ 
tor  PSD  (marked  by  x).  Right  panel,  corresponding  autocorrelation 
functions.  L^  =  1.75  Km,  S  =  -5/3,  a2  =  1.02E-03,  spacing  =  100m. 
Modified  theoretical  autocorrelation  function  at  lag  =  0.  One  predic¬ 
tor  coefficient. 
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12.  Left  panel,  theoretical  PSD  input  (unmarked),  autoregressive  predic¬ 
tor  PSD  (marked  by  x),  and  simulated  PSD  (marked  by  small  square). 
Right  panel,  corresponding  autocorrelation  functions.  Lq  =  10  Km,  S 
=  -5/3,  a2  =  4.9E-03,  spacing  =  100m.  Theoretical  autocorrelation 
function  modified  at  lag  =  0.  Six  predictor  coefficients. 

13.  Autoregressive  simulated  structure  sequence.  Lq  =  10  Km,  S  =  -5/3, 
a2  =  4.9E-03. 

14.  Conclusions.  Self-explanatory. 

15.  The  auto-regression  Yule-Walker  Matrix  equation  for  solving  the  AR 
coefficients. 

16.  The  Levinson  Recursion  Algorithm.  An  iterative  solution  of  the  Yule- 
Walker  equations. 

17.  A  forward  and  backwards  method  for  minimizing  the  error  of  an  AR 
process  power  spectrum. 
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nislograiii  of  sample  Gaussian  random  number  sequence.  Theoretical  mean  =  0..  theoretical  S.D.  =  0.055. 
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Figure^  Left  panel.  Iheoretlcal  F’SD  Inpul  (unmarked),  autoregressive  predlelor  PSD  (marked  by  x).  and  simulated  PSD  (marked 
Q>  by  small  square).  Right  panel,  corresponding  autocoi relation  functions.  L,,  =  1.75  km.  S  =  -5/3.  a  =  1.02E  03. 

'  spacing  =  100  m.  Theoretical  autocorrelation  function  modified  at  lag  =  0.  Six  predictor  coefficients 
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Flgure;^^^  Histogram  of  autoregressive  simulated  structure  sequence.  =  1.75  km.  S  =  -5/3,  o  =  1.02E-03 
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riyure^  l^ft  panel,  theoretical  PSD  Input  (uninarked).  autoregressive  predictor  PSD  (marked  by  x).  Right  panel,  corresponding 
//  autocorrelation  functions.  4  =  1.75  km.  S  =  -5/3  a  =  1.02E  03.  spacing  =  100  m.  Modified  theoretical 
autocorrelation  function  at  lag  =  0.  One  predictor  coefficient 
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Figure  12.  Left  panel,  llieoietical  FSD  Input  (nninarkeci).  autoregressive  predictor  PSD  (marked  by  x).  and  |lmutaled  PSD  (marked 
by  small  square).  Right  pane!,  corresponding  autocorrelation  functions.  =  10  km,  S  =  -5/3,  a  =  4.9E-03,  spacing  = 
100  m.  Tlieorellcal  autocorrelation  function  modified  at  lag  =  0.  Six  predictor  coelTlclenls 
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Figure  13.  Auloreyressive  sliiuilalecl  sluictuie  sequence.  L,.  -  10  km.  S  -  -5/3,  a  4.9E  03 
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Yule  Walker  Equations 
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2(M-NN) 


AN  UPDATE  ON  THE  AFGL  OPTICAL  TURBULENCE 
RADIOSONDE  MODEL 


Edmond  M.  Dewan 


Phillips  Laboratory/  Geophysics 

29  Randolph  Road 

Hanscom  AFB,  MA  01731-3010 


Recently  the  AFGL  model  has  been  applied  by  Lt  Col  Roadcap  of  Phillips  Lab  (at  Kin) and 
AFB)  to  radiosonde  data  for  the  purpose  of  estimating  the  effects  of  turbulence  on  an  AF 
Airborne  Laser  Weapon  System  (ABL).  This  had  the  effect  of  bringing  new  attention  to  this 
model  developed  some  years  ago  by  Dewan,  Good,  Beland,  and  Brown*.  While  it  is  true  that 
our  model  has  been  described  in  the  open  literature,  these  descripdons  left  out  of  account  certain 
important  items  of  interest.  The  purpose  of  this  talk  will  be  to  make  these  available  prior  to  their 
publication.  (An  in-house  report  is  in  press).  These  include  (a)  the  basis  of  the  model  and  how 
it  was  constracted,  (b)  the  estimated  uncertainty  of  its  predictions,  and  (c)  potential  pitfalls  in  its 
use.  An  alternate  (NOAA)  model  exists  and  comparisons  will  be  made  between  it  and  the  AFGL 
Model.  Briefly,  in  contrast  to  the  NOAA  model,  the  AFGL  model  is  (a)  an  order  of  magnitude 
simpler  (i.e.  faster),  (b)  has  no  "adjustable  parameters"  (site  location  dependence),  and  (c)  it  may 
possibly  be  more  reliable  in  its  application  to  the  ABL  program. 

It  has  been  indicated,  with  the  help  of  this  model,  that  jet  streams  may  play  an  important  role  in 
how  any  future  ABL  system  might  be  operated. 


*In  house  report  in  press  "A  Model  for  (Optical  Turbulence)  Profiles  Using  Radiosonde 
Data". 
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DEWAN,  GOOD.  GROSSBARD,  AND  BROWN 


ROLE  OF  IN  OPTICAL  TURBULENCE 

effects  on  laser  beam  propagation 
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ARE  POSSIBLE  FROM  EXISTING  PUBLISHED  DATA. 
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A  NIGHTTIME  STRUCTURE  MODEL  OF  ATMOSPHERIC  OPTICAL 
TURBULENCE,  C^"  DERIVED  FROM  THERMOSONDE  AND 
HIGH  RESOLUTION  RAWINSONDE  MEASUREMENTS 


James  H.  Brown 


Phillips  Laboratory/  Geophysics 

29  Randolph  Road 

Hanscom  AFB,  MA  01731-3010 


Data  from  fifteen  thermosonde  flights  was  used  to  develop  a  simple  nighttime  structure  model 
of  C^.  High  resolution  rawinsondes  provide  fine  scale  estimates  of  atmospheric  temperature 
gradients  and  variances  of  wind  speed.  Non-linear  regression  applied  between  the  thermosonde 
and  rawinsonde  measurements  provide  the  model  estimator.  This  quasi-empirical  model  is 
based  upon  the  theoretical  description  given  by  Tatarski  and  upon  an  exponentially  scaled 
estimate  of  a  theoretical  model  of  the  eddy  dissipation  rate  given  by  Weinstock.  A  discussion 
of  previous  models  and  a  comparison  with  the  Dewan  et.al.  microshear  model  is  presented. 
Model  profiles  computed  for  other  sites  and  seasons  is  compared  favorably  to  related 
thermosonde  profiles. 
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Figure  Captions 

1.  Turbulence  Structure  Constant  equation.  depends  on  temperature 
fluctuations  C{  . 

2.  Cn  Hufnagel  model.  expressed  as  a  function  of  altitude  and  RMS 
wind  speed  only. 

3.  NOAA  Cn  model.  expressed  as  function  of  statistical  distributions 
of  Brunt- Vaisalla  frequency,  wind  shear,  and  outer  scale. 

4.  The  Dewan  et.al.  model.  expressed  as  function  of  outer  scale  where 
outer  scale  is  modeled  after  high  resolution  wind  shears. 

5.  Thermosonde  model.  expressed  as  function  of  Brundt- Vaisalla 
frequency  and  outer  scale.  Outer  scale  modeled  after  thermosonde 
measurements . 

6.  Thermosonde  Model  1  and  Model  2  expressions.  Outer  scale  expressed  as 
functions  of  RMS  wind  speed  and  Brundt-Vaisalla  frequency. 

7.  Final  Form  of  expression  for  outer  scale  model.  Ratio  of  scaled  wind  speed 
variance  and  Brundt-Vaisalla  frequency. 

8.  Thermosonde  Profiles  measured  at  Pennsylvania  State  University,  Flight 
L4007,  May  4,  1986,  Temperature,  Relative  Humidity,  C^,  Wind  Speed, 
Wind  Direction. 

9.  Thermosonde  Cf  measurement  (raw  and  smoothed)  profiles  compared  to 
model  (1)  Cf  profile.  Flight  L4007. 

1 0.  Binned  scatter  plot  of  data  from  entire  Pennsylvania  State  University  cam¬ 
paign.  "L"  for  smoothed  thermosonde  measurements  compared  with  "L" 
from  model  (1).  A  45  degree  slope  represents  perfect  agreement.  Error  bars 
represent  the  standard  deviation  of  the  data  in  each  bin.  Numbers  above 
the  plots  are  the  number  of  points  falling  outside  three  standard  deviations. 
Left-hand  plot  is  troposphere.  Right-hand  plot  is  stratosphere. 

1 1 .  Regression  constants  for  Model  (1)  and  Model  (2). 
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12.  Thermosonde  Profiles  measured  at  Champaign,  Illinois  Flight  LI  01 4,  June 
1988. 

13.  Thermosonde  Cf  measurement  (raw  and  smoothed)  profiles  compared  to 
model  (1)  Cf  profile.  Flight  LI  014 

14.  Outer  scale  -  smoothed  thermosonde  data  compared  to  model. 

15.  Smoothed  shear  profile  as  determined  by  Dewan's  application. 

1 6.  Thermosonde  Cf  measurement  (raw  and  smoothed)  profiles  compared  to 
Dewan,  et.al.  Cf  model  profile.  Flight  L4007. 

17.  L(z)  from  smoothed  thermosonde  data  compared  to  Dewan,  et.al.  L(z) 
model  profile.  Flight  L4007. 

18.  Brunt-Vaisala  frequency  and  RMS  wind  speed  profiles  derived  from 
smoothed  measurement  compared  with  Dewan  et.al.  model  Cf  profile. 
Flight  L4007. 

1 9.  Scatter  plot  of  "L"  for  smoothed  thermosonde  measurements  compared  with 
"L"  from  Dewan,  et.al.  model.  Flight  L4007.  Leftmost  plot  for  troposphere. 
Rightmost  plot  for  stratosphere. 

20.  Isopianatic  Angle  expressed  as  function  of  Cf  and  altitude. 

21.  Isopianatic  angle  calculations  for  Pennsylvania  State  University  Campaign. 
Comparison  of  measurement  to  model. 

22.  Dewan  et.  al.  model.  Isopianatic  angle  calculations.  Comparison  of  meas¬ 
urement  to  model. 

23.  Isopianatic  angle  calculations  for  Pennsylvania  State  University  Campaign, 
Champaign,  III  Campaign,  and  Desert  Site  Campaigns.  Comparison  of 
measurements  to  model. 

24.  Conclusions.  Self-d^criptive. 
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Figure^  'Iliermosoiicle  C,  Measurement  (Raw  and  Smoothed)  Profiles  Compared  to  Motlcl  (I)  C,  Profile.  Fliglil  L4007 
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Pignrti^  Binned  Scalier  Plot  of  Data  from  Rnlirc  Pennsylvania  Slate  University  Cain|)aign.  "L"  for  smoollicd  lliermosonde  measnreinenis  compared  willi 
/O,  “13'  from  Model  ( I ).  A  45  degree  slope  represenis  perfect  agreement.  Error  bars  rc()rescnt  the  standard  deviation  of  llie  data  in  each  bin,  Nnmlrcrs 

above  tlie  plots  are  the  nnmbcrof  |X)inls  falling  outside  three  standard  deviations.  Ixft-liand  plot  is  lroiK)spliere.  Riglit-liand  plot  is  stratosphere 
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Figure  16.  Thermosonde  cf  measurement  (raw  and  Smoothed)  Profiles  Compared  to  Dewan,  et  al.  C,^  Model  Profile.  Flight  L4007 


SiiuK)thccl  niciiiiosoiulc  D;il:«  Compareti  lo  Dewaii,  c,t  at.  l.(/.)  Model  Profile.  Flight  l>l(K)7 
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Measurement  Vs.  Model 
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Dewan  Std.  Dev.  1 .97 


Conclusions 


D) 


0 


E 

CO 

*0 

c 

0 

£ 

Q. 

O 

0 

> 

0 

■D 

0 

0  0 

So 

?  E 

O  (3 

o  .2 

El 
®  o 
^  0 
0 


0 

O) 

0 

C 

0 

> 

■o 

0 

0 

Js£ 

0 


0 

O 

0 

0 


O 

0 


0 


O  _ 

0  CD 

-o 

_  o 

O)  ^ 

ic  o 


T3 

_0 

0 

O 

0 

0 


0 

O 


0 

0 

TJ 


Q  O 

E  B 

0  C 


0 


0 


Q.E 

0 


^0 
m  CO 
ili  0 

V-  Cl 
0  ^ 

3  ® 

C 

15  o 

0  0  0 

’i  E 

|CC  CO  ^ 


0 


0 


O) 


0 

0 

0 


0 

T3 

O 


< 

< 

O 


■ti  O 
^  ■■ 


^  O 
~  0 
~  0 

CO  ^ 

Q-0 

O  O) 

O  Z 


Fidelity  to  structure 


Wednesday  9  June  1993  p.m. 


SESSION  F:  CLIMATOLOGIES 
Chair:  James  H.  Chetwynd,  PL/GPOS 


557 


THE  MOSART  GLOBAL  CLIMATOLOGICAL  AND 
TERRAIN  DATA  BASES 


Dr.  William  M.  Cornette 


Photon  Research  Associates,  Inc. 

10350  N.  Torrey  Pines  Road,  Suite  300 
La  JoUa,  CaLLfomia  92037-1020 


The  MOSART  program  uses  a  number  of  global  climatological  and  terrain  data  bases.  These 
data  bases  include  such  parameters  as  terrain  altitude,  terrain  scene  type,  surface  air  temperature, 
sea  surface  temperature,  snow  cover,  amount  of  terrain,  cloud  cover  (low,  mid,  and  high  etage, 
plus  cirrus).  Attnospheric  data  bases  include  the  six  (6)  MODTRAN  2  model  atmospheres,  plus 
the  additional  seventeen  (17)  APART  model  atmospheres.  These  model  atmospheres  cover  the 
Northern  hemisphere  (equator  to  pole  in  15  degree  increments),  and  include  the  U.S.  and  Israeli 
Standard  atmospheres.  These  data  bases  are  used  to  construct  both  background  definition  and 
atmospheric  conditions  for  any  location  on  the  globe.  This  presentation  will  present  the  data 
bases  and  will  discuss  how  they  are  used  in  the  MOSART  program.  Potential  growth  areas  for 
the  data  bases  and  their  applications  in  the  MOSART  program  will  also  be  presented. 
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Photon  Research  Associates.  Inc. 
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OCEAN  SURFACE  TEMPERATURE 


Ref.:  S.  Leviliis,  Climatological  Atlas  of  the  World  Ocean 
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SHARC/SAMM/MODTRAN  CALCULATIONS  USING  A 
CLIMATOLOGY  MODEL  ATMOSPHERE  GENERATOR 


S=  Adler-Golden,  J.  Gruninger, 
and  M.  Matthew 


Spectral  Sciences  Inc. 
99  S.  Bedford  St. 
Burlington,  MA  01803 


Upper  atmospheric  IR  radiances  have  been  simulated  using  new  input  profiles  that  describe 
ladtude,  seasonal,  solar/geomagnetic,  and  diurnal  (including  solar  terminator)  variabilities.  The 
profiles,  from  the  SHARC  Atmosphere  Generator,  are  based  on  a  combination  of  semi-empirical 
models,  including  the  new  NRL  climatology  database  and  the  MSISE-90  model,  and  theoretical 
calculations.  Agreement  with  field  experiments  is  gready  improved,  and  the  terminator  behaviors 
of  O3  and  OH  emissions  can  be  simulated  with  SHARC  for  the  first  time.  The  profiles  can  be 
generated  down  to  sea  level  for  use  with  SAMM,  MODTRAN,  or  other  codes  that  cover  the 
lower  atmosphere. 
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SAMM/MODTRAN  Trace  Species 


SHARC/SAMM  ATMOSPHERE  GENERATOR 
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Output  Formats  Compatible  with 
MODTRAN/LOWTRAN 


SHARC/SAMM  ATMOSPHERE  GENERATOR 
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IR  RADIANCE  CALCULATIONS 
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TERMINATOR  BEHAVIOR 
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Oo  TERMINATOR:  DATA/MODEL  COMPARISON 


SZA  (degrees) 
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DEVELOPMENT  OF  A  SIGNAL-TO-NOISE  PERFORMANCE  MODEL 

wrrrnN  backscat  version  4.o 


M,G.  Cheifetz:,  D.R.  Longtin,  and  J.R.  Hummel 


SPARTA,  Inc. 

24  HartweU  Ave. 
Lexington,  MA  02173 


The  lidar  backscatter  simulation,  BACKSCAT  Version  4.0,  has  expanded  its  capabilities  and 
now  includes  a  comprehensive  and  versatile  signal-to-noise  performance  model.  This  signal-to- 
noise  (SNR)  model  can  give  performance  and  sensitivity  estimates  for  both  direct  detection  and 
coherent  (heterodyne)  lidar  systems.  The  model  contains  all  the  important  noise  sources  inherent 
in  the  detection  process  and  allows  various  detector  types  to  be  simulated  and  analyzed.  In  this 
paper  we  will  summarize  the  SNR  model  and  its  inputs,describe  some  of  the  buUt-in  detectors, 
and  give  examples. 

Research  Supported  by  Phillips  Laboratory,  Geophysics  Directorate 
Contract  F19628-91-C-0093 
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DEVELOPMENT  OF  A 
SIGNAL-TO-NOISE  PERFORMANCE 
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Work  Performed  Under  Contract  F19628-C-91-0093 
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BACKSCAT  VERSION  4.0 
A  LIDAR  SIMULATION  TOOL 
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BACKSCAT  SIGNAL-TO-NOISE 
PERFORMANCE  MODEL 
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SNR  PERFORMANCE  MODEL 
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Receiver  Field-of-View 


SNR  MODEL  ASSUMPTIONS 
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MODEL  DETECTORS 
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MODEL  DETECTOR 
RESPONSITIVITES 
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ATMOSPHERIC  PARAMETERS 
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PROFILE 


PERFORMANCE  PREDICTIONS 
-  DIRECT  DETECTION  - 


Range/Altitude  (km) 
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AN  OPTICAL  PROFILE  FUNCTION  FOR  MODELING  EXTINCTION  AND 
BACKSCATTER  COEFHCIENTS  IN  VERY  LOW  STRATUS  CLOUDS 

AND  SUBCLOUD  REGIONS 


Neal  Kilmer 


Henry  Rachele 


Physical  Science  Lab 
New  Mexico  State  Univ. 
Las  Cruces,  NM  88003 


Battlefield  Envir’ment  Directorate 
U.S.  Army  Research  Lab 
WSMR,  NM  88002 


A  theoretically  based  microphysics  model  developed  by  the  authors  simulated  135  vertical 
profiles  of  drop  size  distributions  in  and  below  very  low  stratus  clouds.  These  profiles  were 
computed  as  a  function  of  air  mass  type,  maximum  liquid  water  content,  and  surface  (2  m) 
values  of  temperature,  relative  humidity,  and  visibility  representative  of  worldwide  conditions. 
These  drop  size  distribution  profiles  with  Mie  efficiency  factors  simulated  vertical  profiles  of 
extinction  and  backscatter  coefficients  for  eight  wavelengths.  The  extinction  and  backscatter 
profiles  were  fit  with  the  Rachele-Kilmer  (RK)  optical  profile  function  to  significantly  simplify 
computation  procedures.  All  constants  required  to  evaluate  the  RK  optical  profile  function  have 
been  placed  in  computer-accessible  storage,  and  a  computer  program  for  performing  the 
calculations  has  been  prepared  and  is  available  to  DoD  users.  Use  of  this  program  is  described. 
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SOME  PRODUCTS  OF  RK  MODEL 


Liquid  water  content 
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EXnNCnON  COEFFICIENT, 
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SJ0j0UJ 


PROGRAM  RKOPF  -  USER  OPTIONS 
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Type  of  profile  (extinction  and/or 
backscatter  coefficient) 


FORM  OF  OUTPUT  FILE 
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above  ground  level  (both  real 
using  free  format) 
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profiles 


CONCLUSIONS 
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(192.67.8.5)  using  telnet,  logging  on  as  rkopf, 
and  then  entering  rkopf  (lower-case  letters). 
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TIME  AND  POLARIZATION  DEPENDENT  DOUBLE  SCATTERING 
CALCULATIONS  OF  LIDAR  RETURNS  FROM  WATER  CLOUDS 


Richard  Garner 


PhotoMetrics,  Inc. 
Wobum,  MA  01801 


We  have  developed  a  double  scattering  lidar  program  which  is  used  to  calculate  lidar  returns 
from  water  clouds.  The  program,  which  is  implemented  on  a  PC,  is  used  to  interpret  data 
acquired  with  the  Air  Force  Phillips  Laboratoiy’s  (Geophysics  Directorate/GPOA)  elastic 
backscatter,  polarization  sensitive,  Nd:YAG  based  lidar  system.  The  program  determines  the  four 
Stokes  parameters  of  the  backscatter  lidar  radiation,  from  spatially  inhomogeneous  media 
composed  of  spherical  particles,  as  a  function  of  time  and  as  a  function  of  telescope  focal  plane 
location.  We  use  the  program  to  determine  particle  size  distributions  and  multiple  scatter 
corrected  extinction  coefficients  of  water  clouds.  In  this  talk  we  will  describe  the  program, 
compare  its  results  to  lidar  data,  and  present  and  show  examples  of  our  data  analysis  techniques. 

Work  Supported  by  the  AP  Phillips  Laboratory',  Geophysics  Directorate,  Hanscom  APB,  MA 
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Tiine5  Polarization,  and  Focal  Plane  Dependent 
Calculations  of  Doubly  Scattered  Lidar  Rettirn 

Radiation  from  Inbomogeneons  Water  Clouds, 


Richard  Garner 
PhotoMetrics.  Inc. 

4  Arrow  Dr. 
Woburn,  MA  01801 


Work  supported  by 

U.S.  Ak  Force  Philips  Laborato-ry/ Geophysics  Directorate/GPOA 

Hanscom  AFB,  MA. 


Animal  Reriew  Conference  on  Atstos^feL-eric  Tran^ission  Mo<Jels 
Phill4>3  Laboratory,  Hanscom  ATB,  Massack-isetts 
8-9  June  1993 
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Using 


•  a  doubled  NdiYAG  elastic  backscatter  lidar  (witk 
capability  for  simultaneous  polarizatio;n  me^^mem  and 

•  a  time,  polarization,  and  focal  plane  dependent  double 
(Mie)  scartering  model  (implemented  on  a  PC) 

Do  the  following: 

•  Investigate  regimes  of  validity  of  model, 

•  determine  multiple  scatter  corrections  to  li-dar  derived 
extinction  coeScient,  and 

•  determine  particle  size  distributions. 

Outline  of  this  talk: 

•  Describe  model. 

•  Compare  azimuthally  dependent  lidar  retums  to 
model  calculations. 

•  Describe  and  show  example  of  teckniqlue  to  detenffiEte  size 
distribution  and  multiple  scatter  corrected  extisetio-m 
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Differences  with  Previous  Models 

®  Variable  extinction  along  z  direction 
(e,g.,  extinction  calculated  from'Hdar  data) 

®  Lidar  return  versus  focal  plane  location 
(azimuth  and  polar  angle). 

®  General  incident  Stokes  vector. 

®  Receiver  response  a  function  of  polar  angle. 

®  Nomntegrable  singularity  not  encountered 
(different  variables  of  integration). 

®  Laser  has  zero  divergence. 
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Received  power  from  one  set  of  scatters: 

Z[-6)P{e2)hh)h-'p)Po 

rotation  matrix  (in  space  of  Sto^iES  Tector) 

scattermg  picase  matrix 
®  extrsctic®  coefficient  at  scatter 
®  r:  total  optical  patlli  iengtb 

®  AriO):  receiver  response  (normalized  to  aperture  area) 


Intesration  of  dP, 


•  dPr  is  integrated  over  all  paths  of  the  same  length 
(qniadxuple  integral). 


extinction  coefficient  <7  is  a  function  of  z, 


is  the  Mie  scatter  phase  matrix: 


P{9)^ 


A[e)  0  0  0 

0  B[9)  0  0 

0  0  C{e)  -D[B) 

0  0  D{B)  C{B) 

\ 

(size  distribution  independent  of  position) 
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The  field  of  view  is  a  section  of  an  annular  region 


detection  region 


Integration  of  dPr  [cont,) 


•  Tke  6  integral  can  be  performed  analytically. 

•  The  remaming  ihiee  integrals  (in  Ri.  9i,  R2)  are  performed 
numericallv. 


•  The  general  expression  is  long  and  is  not  presented  here. 


•  When  the  field  of  view  is  a  complete  annular  region 
(0  <  o  <  27r): 


^  ^  /•ct/2  [hhW  Q  Jn  ^ 


^{AiA2  4-  BiB2){^P\\o  d"  B±o)  —  y(C'iC2  —  DiD2){P\\o  ~  P J_o) 
j{AiA2  -T  BiB2){P\\o  "T  ^P±o)  4-  f(C'iC2  —  DiD2){P\\o  —  P±o) 

7rUo{CiOz  —  D1D2) 

27rV,{CiC2  -  D1D2) 
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®  extiactian  as  a  fimction  of  range 
@  particle  size  distribution  (independent  of  range) 
@  complex  index  of  refraction 


®  incident  Stokes  vector 


®  field  of  view  (polar  and  azimutbal  angle  ranges) 
®  aperture  area  and  response  function 
®  temporal  pulse  width 
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[ethod  of  Calculation: 


•  Mie  scatteriag 

—  Calculate  pkase  matrix  coefficients  at  discrete  angles. 
—  Store  in  arrays. 

—  Determine  values  ax  other  angles,  as  needed,  by 
interpolation. 


•  Double  scattering 

—  Calculate  triple  integral  using  three  nested  Simpson's 
rule. 


Calculations  are  imp-lemented  o-n  a  PC. 
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Base  case: 

•  extinction:  10  'km~^  (constant) 

•  rx  =  1.5  m 

•  modified  gamma  size  distribution  with 

Oo 

—  a  —  2 

•  outer  boundary  of  field  of  view:  10  mrad. 

(full  angle) 

•  inner  boundary  of  field  of  view:  0  mrad. 

•  index  of  refraction:  1.33-ril0“® 

•  cloud  base:  1  km 

All  results  shown  are  for  variations  of  one  or  two 
parameters  from  the  base  case. 
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Do.ubie  Scoitering  Mode! 
moafied  gamma;  (Q=?2,xm=35) 
extinction=10  1/km 


scan 


ms 


2  =  0.385 


Cumuius  douc  size  dstribution  mode! 
(Hcncb.ccK  of  Gecphysics  and 
the  Space  Environment,  p.  18—29) 
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Mie  scattering  phase  matrix 
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Multiple  BCAttat*  '6f 

;  Lidar  returns  fi'bill  6<(fpMiUs 

Double  M.o  aoaUer  program  reaulls  (averagud  UVbr  ■l-d:2  'km) 

VS.  azimuth  aittfls 
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Determiniiig  size  distribution  aad 
multiple  scatter  corrected  extinetian 


1.  Acquire  lidar  data  (parallel  and  perpendicular  polarization 
components  of  backscattered  radiation  from  water  cloud) . 
Calculate  depolarizaiion  versus  range  in  cloud. 


2.  Relatively  calibrate  the  two  channels  using  background  daia 
(assumed  to  have  depolarization  of  unity). 

3.  Calculate  the  extinction  coefficient  as  a  function  of  range 
with  the  Klett  technique  using  the  parallel  polarization 
component. 


4.  Use  the  calculated  extinction  profile  in  the  double  scatter 
mcKiel.  Run  the  model  many  times  for  a  variety  of  size 
distributions.  Assume  a  particular  size  distribution 
(e.g.j  modified  gamma)  with  free  parameters  that  can  be 
varied. 

5.  Choose  the  best  size  distribution  by  matching  predicted 
depolarization  to  measured  depolarization. 

6.  Calculate  co'rrection  to  Klett  extinction  above  using  model 
results. 


7.  Go  to  step  4. 
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•  prkn-ary  and  dpu-bled  Nd:YAG 

•  25  mJ/puis€  @  532  nin 

•  20  Hz  pulse  repetition  rate 

•  15  cm  aperture  diameter 

•  10  mrad  field  of  view 

•  20  MHz  (7.5  m)  maximum  da,ta  rate 

•  steerable  (upper  fiemisphere) 

•  two  simultaneous  linear  polarization 
measurements  (532  nm  only) 

(parallel  and  perpendicular  to  laser  polarization) 


'PL/GD/GPOA  Transportable  Optical  Atmospheric  Data  System  (TOADS) 
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powej:  (qrb.  mits) 


u 

n: 

b 

■5 

CL 

0) 

-u 


model/data  comparison 


Q  I  I  I  1  I  t  t  I  t  r  I _ I _ [ _ I _ ■  ■  «  1  « _ I — t — I  I  I  I  I  1  J 

6£  6£5  6.7  6.75  6.S  635  6i9 


range  (km) 
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T?echjg.j.^He  to  determine  multiple  scatter 
carrec.ti:on  to  ext-m-ct-ian: 

Utilize  ideas  presemed  by 

S.R.  Pal  and  A.I.  Carswell. 

/■ 

Apvl  Opt  8,  1990-1995,  (1976) 


Say 

Then 


•^ii^  =  ;^.^exp{-2r 
A 

P\\t  =  exp  1-2(7  -  '„) 


P\lm  =  P\\t  —  P\\..  =  O' Pi 


l|i 


P:1 


Pit  —  aPi 


dr, 


m 


m 


dR 


a  CT' 


m 


®  Measure  Pu  and  Pi  with  hdar. 

IP 

•  Determine  a  wixh  model. 

•  Calculate 


704 


raturn/incldant  power 


Double  scatter  mo 


&85 


>  as 
[km) 


6^ 


double  perp/dout^le  para 


modified  aamma  distribution 
xm=75,  a==3 


r- 

h 


r 

L 

y .  I  I  '  !  I  I  I  I  *  '  '  !  1  I  I  1 _ ^ — I _ L _ I _ : - ! - 

6.6  6.65  6.7  5.75  6.8  6.85 

rcnge  (km) 

extinction  frorri  7/7/92,  =2 


6.-9 
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return/incident  power 


Conclusions 


For  constant  extinction  the  depolarization  of  the  received  radiation 
(single  and  double  scattered)  increases  from  zero  to  some  roughly 
steady  state  value.  The  value  is  relatively  low  (typically  in  the  range 
of  0.1  to  0.2),  even  for  large  extinctions  or  large  fields  of  view. 


The  steady  state  depolarization  increases  with  increasing 

—  extinction 

—  field  of  view 

—  cloud  base 


Most  double  scattered  power  is  from  near  the  axis  where  much  of  it  is 
polarized  parallel  to  the  incident  radiation.  This  is  the  reason  for  the 
relatively  low  steady  state  depolarization. 


Relatively  high  steady  state  depolarizations  are  achieved  for  clouds 
that  are  farther  away  from  the  lidar  (because  regions  farther  from  the 
axis  then  contribute  more  to  the  returned  double  scattered  power). 


For  real  clouds  it  is  important  to  consider  the  varying  extinction  as  a 
function  distance  into  the  cloud.  We  see  agreement  between  the  model 
and  lidar  data  from  clouds  using  the  extinction  calculated  from  the 
data  (using  single  scattering  assumptions). 


We  see  agreement  between  model  and  data  (lidar  and  cloud  chamber) 
as  a  function  of  azimuth  angle. 


The  parallel  component  of  the  backscattered  radiation  versus  azimuthal 
angle  is  sensitive  to  size  parameter. 


SKY  RADIANCE  AND  RAY  BENDING  CALCULATION 


Gertrude  H.  Komfeld 


Modeling  and  Simulation  Branch 
Army  Research  Laboratory 


An  interface  of  MODTRAN  results  with  computer  generation  of  thermal  imagery  is  planned. 
The  misuse  of  MODTRAN  information  might  lead  to  dangerous  conclusions. 

A  closed  from  solution  for  altitude  dependent  sky  radiance  that  considers  the  curvature  of  the 
earth  and  ray  bending,  was  originated  in-house  by  the  SSI  Branch  at  ARL,  and  tested  for  the 
special  case  of  a  Boltzmann  distribution  of  particle  densities  in  the  atmosphere. 

Ray  bending  is  very  important  because  it  potentially  invalidates  laser  designators  and  also  is 
an  indication  of  unstable  potentially  dangerous  atmospheres.  A  flat  earth  approximation  would 
imply  erroneous  total  reflections.  The  method  can  be  amended  for  any  analytically  described 
atmosphere. 
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Solution  by  Law  of  Cosines 

Layering  with  Bent  Rays 

Broad  Band  Transmission  Calculation 
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Total  Reflection  occurs  at  some 
Layer,  in  case  atmosphere  is 
very  unstable 
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imply  a  sine  that  exceeds  unity 
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Corners  at  Earth  Center,  Earth  Surface  on  H 


Closed  Form  Solution 
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Sky  Radiance 
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^UT  »  radiant  powers  at  (l^  ,ly  )  tocationjspecfal  case  the 

sky  has  (1^  ,ly)  according  to  3°K 
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THERIVIAL  ENERGY 
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Prepare  Standard  Atmospheres 

Ray  Bending 


INTEGRATION  OF  LOWTRAN  INTO  GLOBAL 
CIRCULATION  MODELS  FOR  OBSERVING  SYSTEM  SIMULATION 

EXPERIMENTS 


S.A.  Wood  and  G.D.  Emmitt 


Simpson  Weather  Associates,  Inc. 
809  E.  Jefferson  St. 
Charlottesville,  VA  22902 


The  LAWS  Simulation  Model  (LSM)  simulates  observations  from  a  space-based  Doppler  lidar 
wind  sounder.  A  main  component  of  LSM  is  its  atmosphere  generator  model  that  produces 
global  estimates  of  aerosol  optical  properties,  opaque  clouds  and  subgrid  scale  turbulence  using 
output  from  the  European  Center  Medium  Range  Weather  Forecast  (ECMWF)  global  circulation 
model.  A  major  issue  that  will  be  discussed  is  the  reasonableness  of  the  /3  backscatter  fields 
resulting  from  the  integration  of  LOWTRAN  into  the  Global  Circulation  Models  (GCMs). 
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INTEGRATION  OF  LOWTRAN  INTO 
GLOBAL  CIRCULATION  MODELS 
FOR  OBSERVING  SYSTEM 
SIMULATION  EXPERIMENTS 


S.A.  Wood 
G.D.  Emmitt 

Simpson  Weather  Associates,  InCo 
Charlottesville,  VA  22902 


FIGURE  1 
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Figure  1.  Introduction  slide  for  presentation  "Integration  of 
Lowtran  Into  Global  Circulation  Models  for  Observing  System 
Simulation  Experiments"  by  S.  A.  Wood  and  G.  D.  Emmitt. 

Figure  2.  Overview  slide  for  "Integration  of  Lowtran  Into 
Global  Circulation  Models  for  Observing  System  Simulation 
Experiments"  by  S.  A.  Wood  and  G.  D.  Emmitt. 

Figure  3.  Block  diagram  for  the  LAWS  Simulation  Model  - 
Global  Version. 

Figure  4.  Overview  slide  for  the  Observing  System  Simulation 
Experiments  for  LAWS. 

Figure  5.  Overview  slide  for  the  1°  X  1°  Global  LOWTRAN  Input 
Data  Base. 

Figure  6.  Flow  diagram  for  the  LSM's  optical  property  model. 

Figure  7.  LAWS  baseline  signal  to  noise  equation. 

Figure  8.  Global  9.11  jLtm  relative  aerosol  backscatter  (jLtm**2 
m'‘  sr'‘)  at  the  earth's  surface  for  1/16/79  0600Z.  The  aerosol 
backscatter  has  been  multiplied  by  the  lidar  wavelength  squared. 

Figure  9.  Global  2.1  iJ.m  relative  aerosol  backscatter  (jtxm**2 
m'*  sr'*)  at  the  earth's  surface  for  1/16/79  0600Z.  The  aerosol 
backscatter  has  been  multiplied  by  the  lidar  wavelength  squared. 

Figure  10.  Global  relative  humidity  field  at  the  earth's 
surface  for  1/16/79  0600Z. 

Figure  11.  Attenuated  global  9.11  ^i^a  relative  aerosol 
backscatter  (/im**2  m"'  sr'‘)  at  the  earth's  surface  for  1/16/79 
0600Z.  The  aerosol  backscatter  has  been  multiplied  by  the  lidar 
wavelength  squared. 

Figure  12.  Attenuated  global  2.1  jum  relative  aerosol 
backscatter  (/im**2  m'‘  sr'‘)  at  the  earth's  surface  for  1/16/79 
0600Z.  The  aerosol  backscatter  has  been  multiplied  by  the  lidar 
wavelength  squared. 

Figure  13.  Global  9.11  juro  relative  aerosol  backscatter  (jLtm**2 
m'^  sr'‘)  greater  than  2.5  e-6  with  integrated  cloud  cover  less  than 
90%  for  1/16/79  0600Z.  The  aerosol  backscatter  has  been  multiplied 
by  the  lidar  wavelength  squared. 

Figure  14.  Global  integrated  cloud  cover  from  the  top  of  the 
atmosphere  to  the  earth's  surface  for  1/16/79  0600Z. 
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OBSERVATION  ERRORS 


LOG  8ACKSCATTER  (M-l  SR-I)  AT  850  KM 


FHE  OBSEFIVING  SYSTEM  SIMULATION  EXPERIMENTS  FOR  1  AWS 


nos  ECkfWF 
MAfURE  BUN 


LOMHAHt 

WITH 

nature  suh 


OSSE 


(  LAWS  I  i 
iaiMULATIC^ 
l|  MODEL i 


IGLOOAL  VERSIOMI 


LSM 

I  OUTPinS 


1^x1°  GLOBAL  LOWTRAN  INPUT 

DATA.  BASE 


LOCATION  PROFILE 

-  Tropical 

-  Subtropical 

-  Midlatitude 

-  Sub-Artic 

•  HAZE  MODEL 

-  Rural 

-  Navy  Maritime 

-  Ocean 

-  Urban 

-  Tropospheric 

-  Desert 

•  COASTAL  INFLUENCE 

“  Open  Ocean 

-  Weak  Continental  Influence 

-  Strong  Continental  Influence 


FIGURE  5a 
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STRATOSPHERIC  MODEL 

-  Background 

-  Moderate  Aged  Aerosol 

-  Moderate  Fresh  Aerosol 

-  High  Aged  Aerosol 

-  High  Fresh  Aerosol 

Upper  Atmosphere  Model 

-  Normal  Upper  Atmosphere 

-  Extreme  Upper  Atmosphere 

-  Transition  Volcanic  to  Normal 

-  Transition  Norma!  to  Volcanic 


f(  WAVEL,  RH,  TEMP.,  PRESS.,  Z) 


fO/lS(RH),  Z,  SEASON,  PRCFiLE) 


f(RH,  WAVELENGTH) 


f(  WAVEL.,  RH,  WIND,  COASTAL  INF.) 


F  —  GuRE  6s, 
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Baseline  LAWS  SNH  Equation 


SNE^ 


TT 


•ni-n2-n3-n4-J*D"-X-B-e-^^“®‘'" 


8*h  B-R^ 


=  quantum  efficiency  =  .40 

=  optical  efficiency  =  .65 

xij  =  system  efficiency  racuor  =  .32 

n^  =  other  losses  =  .  5 

j  =  laser  energy  (Joules) 

D  =  mirror  diameter  (m) 

X  =  laser  wavelength  (m) 

B  =  backscatter  (m“  sr‘^) 

Q,  =  atuenuation  as  a  luncuion  o;. 
telescope 

h  =  Planck's  constant  =  6.63  x  10 
3  =  mrocessing  bandwidth  =  2  V^/ X 

V5  =  search  window  (ms*)  =  50  =  (- 

R  =  slant  range  (m) 


FIGURE  7 


=  wideband  SNR 
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INTrGRATED  CLOUD  COVER  FROM  THE  TOP  OF 
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AEROSOL  DISTRIBUTION  AND  IR  BROADBAND  TRANSMITTANCE 
IN  THE  MARINE  BOUNDARY  LAYER  IN  THE  MEDITERRANEAN 
ENVIRONMENT 


Mireille  Tanguy,  Michel  Autric,  and  Bernard  Salles 


DCN/CESDA  Mourillon/GR.  OPT 
BP  77 

83800  Toulon  Naval  FRANCE 


An  IR  broadband  tranmissometer  has  been  settled  near  the  Toulon’s  coast  in  the  Mediterranean 
Sea.  The  atmospheric  transmittance  has  been  measured  along  a  horizontal  path  of  8  km  at  about 
30  m  above  the  sea  surface.  Meteorological  parameters,  aerosol  density  distribution  and 
visibility  have  been  recorded  too.  Partial  results  have  been  presented  during  the  SPIE  meeting 
in  Orlando  (April  1991).  So  further  conclusions  will  be  presented  on: 

-  Aerosol  density  behavior 

-  Atmospheric  transmittance  behavior 

-  Systematic  comparison  between  measures  and  LOWTRAN  7 

using  Navy  maritime  model  and  Maritime  model. 
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DCN  TOULON  (TRANCE) 

Directorate  for  Naval  Construction  working  for  the  French  Navy 
CTSN  /  LSA  /  GR.  OPTRONIQUE 
Technical  Center  of  Naval  Systems 


LM.F.  of  MARSEILLE  (FRANCE) 
Fluids  Mechanisms  Institute 


author  :  Miss  Mireille  TANGLTY 


co-authors :  M.  Herve  BONHOMMET  (DCN) 

M.  Bernard  SALLES  (DCN) 

M.  Michel  AUTRIC  (IMF) 
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OUTLINE 


INTRODUCTION 


THE  GEOGR.\PHIC  SITE  AND  THE  TECHNICAL  MEANS 


THE  AEROSOL  DENSITY  MEASUREMENTS 


THE  ATMOSPHERIC  TRANSMITTANCE  MEASUREMENTS  AND 
THE  LOWTRAN  7  CODE 


CONCLUSION 
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J 

INTRODUCTION 


FINAL  GOAL  : 


-  Improvement  of  the  software  :  PREDIR  used  for  the  prediction  of  the 
range  of  the  Navy's  ER  optronics  passive  systems 

-  LOWTRAN  code  :  subroutine  of  PREDIR  for  radiance  and 
transmittance  calculations 


PIER  PRESENT  STUDY: 

-  Validation  of  the  LOWTRAN  version  7  in  a  marine  mediterranean 
environment  for  the  following  conditions ; 

Atmospheric  Transmittance 

aerosol  contribution  at  14  meters  high 
horizontal  path  at  30  meters  high  in  a  coastal  zone  ^ 
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GE0GK4PHIC  SITE  AND  THE  TECHNICAL  MEANS 


iiS: 


opticalcounters 


Yisibility  meter 


emission: 


EAST 


TOULON 


WEST 


meteo  station 


's  roads 


station 


period 


AEROSOLS 


Aout-September 


TRANSMITTANCE 


Oct-Nov-Dec 


CSASP  100  HV  14  m  (building) 

ASASP-X  14  m  (building)  14m(building)  -  7m(ship) 

TRANSMSSOMETER  3-5  /  8- 12pm  -  30m  -  8km 

METEOFRANCE  125  m  (2km  to  Soutii) 

VISIBILITY  METER  17  m  (building)  17m(building)  -  8m  (ship) 

AUTOMATIC  METEO  40m  (receiving  set) 


ASASP-X 


AEROSOL  DENSITY  MEASUREMENTS 


D(um) 


3  kinds  of  exploitation  : 


averaae  volume  over  4  diameters  ranges  ; 

0.09um-0.5um  0.5pm-2um  2pm-10!im  10um-47pm 


volume  d\  —  d2 


(^2 


rdl 


4  3  dN  , 

—  X  TUX  r-^  X  — —  X  dr 

3  dr 


extinction  coefficients  Kaer  calculated  from  measured  aerosols 
distributions  (Mie  theory) 

total  transmittance  by  introducing  extinction  coefficients  into  user 
defined  card"  of  Lowtran  code 
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AEROSOLS  DENSITY  MEASUREMENTS 


The  average  volumes  and  the  meteorological  conditions 


bad  relation  between  air  mass  and  [O.ljim-O.Sum]  aerosol  density 

-  proximity  of  coast 

-  likeness  of  the  trajectory  of  depression  and  anticyclone  air  mass 


volume  0.1-0. Sum 
\  ‘ 


AEROSOLS  DENSITY  MEASUREMENTS 
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The  average  volumes  and  the  meteorological  conditions 


AEROSOL  DENSITY  MEASUREMENTS 
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aerosols  extinction  coefficients 


ATMOSPHERIC  TRANSMITTANCE  MEASUREMENTS 


9 


accuracy  of  temperature  measurements 

feature  of  thermometer '  s  location 

aerosols  density 's  variations  are  not  preponderant 
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ATMOSPHERIC  TRANSMITTANCE  MEASUREMENTS 

LOWTRAN  7  CODE 
NAVY  AEROSOL  MODEL 


black  points 
white  points 


using  visibility 
without  visibility 
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ATMOSPHERIC  TRANSMITTANCE  MEASUREMENTS 


LOWTRAN  7  CODE 
MARITIME  MODEL 


754 
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CONCLUSION 


COASTAL  MEDITERRANEAN  ZONE 


Navy  Aerosol  Model  unadapted  at  14  meters  high 


reference  aerosol  distribution 

surface  range 
+ 

absolute  humidity 

estimation  of  atmospheric  transmittance  at  30  meters  high 
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MARrriME  MODIFICATIONS  TO  LOWTRAN  RADIANCE 


i 


P 

C.R.  Zeisse 


NCCOSC  RDTE  DIV  543 
53170  Woodward  Road 
San  Diego,  CA  92152 


LOWTRAN  6  has  been  modified  for  shipboard  observation  (HI  =  10  m)  of  the  marine  horizon 
(ANGLE  =  90  ±  1°)  in  the  long  wave  band  (830  to  1250  cm'^).  For  paths  to  the  sky  an 
anomalous  dip  originally  occurred  in  the  radiance  calculated  just  above  the  horizon.  This  dip 
disagreed  with  low  altitude  observations  of  the  maritime  sky  and  was  removed  by  increasing  the 
atmospheric  layering.  For  paths  to  the  earth,  the  earth  has  been  reinterpreted  as  the  sea.  Using 
Cox-Munk  wave  slope  statistics  for  the  sea  surface,  the  following  sea  radiance  contributions  have 
been  introduced  in  addition  to  the  path  radiance  already  provided  by  LOWTRAN  6;  (1)  thermal 
emission  from  the  sea,  (2)  reflection  of  sky  radiance  by  the  sea,  and  (3)  solar  glints.  These 
modifications  increase  the  calculated  radiance  by  as  much  as  a  factor  of  two  and  bring 
calculations  to  within  a  few  °C  of  marine  observations.  They  are  being  considered  for 
introduction  into  future  LOWTRAN  versions.  Finally,  it  will  be  proven  that  the  Ben-Shalom 
radiance  formula  does  not  respond  to  aerosol  content,  making  it  inappropriate  for  use  in  the 
marine  environment. 
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Advanced  Navy  Aerosol  Model 
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9  Dec  1986 


ZENITH  ANGLE  (degrees) 

Comparison  of  measured  sky  radiances  and  those  calculated 
by  LOWTRAN  6  vs  zenith  angle. 


RIZONTAL  GEOMETRY  WITHIN  LAYER  1 


30  June  1987 
8-12  [im 
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Fig.  4.  Comparison  of  sky  radiances  measured  with  the  AGA  sys¬ 
tem  on  30  June  1987  with  radiances  calculated  with  LOWTRAN  6  with 

and  without  layers  added. 
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OUTLINE  OF  COX-MUNK  APPROACH 
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VERTICAL  RADIANCE  SLICE;  10  OCTOBER  1991 


Radiance  (mW  str 


0  OCTOBER  1991  HORIZON  RADIANCE 
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10OCTV2 


10  OCTOBER  1991  HORIZON  RADIANCE  (b  to  12  iim) 
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LOWTRAN  6  Radiance  (Lowtran  4  manual,  p.  9) 
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N*iT) 


For  horizontal  paths  Ben-Shalom  gives  N*(T),  the  blackbody 


OCTOBE 
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VERTICAL  RADIANCE  SLICE:  10  OCTOBER  1991 
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Radiance  (mW  cm'^  str'^) 


.  Maritime  modifications  to  LOWTRAN  6  improve 
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inappropriate  for  marine  use. 
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